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1 INTRODUCTION 


In high energy collisions, perturbative Quantum Chromodynamics (pQCD) successfully predicts inclu¬ 
sive energy spectra of particles in jets. They have been determined within the Modified Leading Log¬ 
arithmic Approximation (MLLA) d O as functions of the logarithm of the energy (ln(l/.xj) and the 
result is in nice agreement with the data of - e + e~ and hadronic - colliders and of deep inelastic scat¬ 
tering (DIS) (see for example 0 0] (5| ). Though theoretical predictions have been derived for small x 
(energy fraction of one parton inside the jet, x 1) 4 , the agreement turns out to hold even for x ~ 1. 
The shape of the inclusive spectrum can even be successfully described by setting the infrared transverse 
momentum cutoff Qq as low as the intrinsic QCD scale A qcd (this is the so-called “limiting spectrum”). 

This work concerns the production of two hadrons inside a high energy jet (quark or gluon); they 
hadronize out of two partons at the end of a cascading process that we calculate in pQCD; considering 
this transition as a “soft” process is the essence of the “Local Parton Hadron Duality” (LPHD) hypothesis 
m m m. that experimental data have, up to now, not put in jeopardy. 

More specifically, we study, in the MLLA scheme of resummation, the double differential inclusive 1 - 
particle distribution and the inclusive k± distribution as functions of the transverse momentum of the 
emitted hadrons; they have up to now only been investigated in DLA (Double Logarithmic Approxima¬ 
tion) £Q. After giving general expressions valid at all x, we are concerned in the rest of the paper with 
the small x region (the range of which is extensively discussed) where explicit analytical formulae can be 
obtained; we furthermore consider the limit Qq ~ A qcd, which leads to tractable results. We deal with 
jets of small aperture; as far as hadronic colliders are concerned, this has in particular the advantage to 
avoid interferences between ingoing and outgoing states. 

The paper is organized as follows: 

• The description of the process, the notations and conventions are presented in section 0 We set there 
the general formula of the inclusive 2-particle differential cross section for the production of two hadrons 
h\ and /y at angle 0 within a jet of opening angle ©o, carrying respectively the fractions x\ and jy of 
the jet energy E\ the axis of the jet is identified with the direction of the energy flow. 

• In sectional we determine the double differential inclusive 1-particle distribution t n n e * C)I 

the hadron h\ emitted with the energy fraction x\ of the jet energy E, at an angle 0 with respect to the 
jet axis. This expression is valid for all x; it however only simplifies for x <C 1, where an analytical 
expression can be obtained; this concerns the rest of the paper. 

• In section 0 we go to the small x region and determine d\n(ifx{) d.\n<s> > X] ^ 1 both for a gluon jet 
and for a quark jet. It is plotted as a function of In k± (or In 0) for different values of £\ = 1 n(1 /a: i); 
the role of the opening angle ©o of the jet is also considered; we compare in particular the MLLA 
calculation with a naive approach, inspired by DLA calculations, in which furthermore the evolution of 
the starting jet from ©o, its initial aperture, to the angle 0 between the two outgoing hadrons is not taken 
into account. 

The MLLA expressions of the average gluon and quark color currents < C > g and < C > q involve 
potentially large corrections with respect to their - expressions at leading order; the larger the (small) x 
domain extends, the larger they are; keeping then under control sets the bound t = In ^ > 2.5. 

4 as the exact solution of the MLLA evolution equations 
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• In section |5] we study the inclusive k± distribution d , which is the integral of dx f d ^ ne with 
respect to xq; It is shown in particular how MLLA corrections ensure its positivity. The domain of 
validity of our predictions is discussed; it is a k± interval, limited by the necessity of staying in the 
perturbative regime and the range of applicability of our small x approximation; it increases with the jet 
hardness. The case of mixed gluon and quark jets is evoked. 

• A conclusion briefly summarizes the results of this work and comments on its extensions under 
preparation. 

Five appendices complete this work; 

• Appendix lAlis dedicated to the MLLA evolution equation for the partonic fragmentation functions 
Dg or q and their exact solutions fUO. They are plotted, together with their derivatives with respect to 
ln(l/re) and In k±. This eases the understanding of the figures in the core of the paper and shows the 
consistency of our calculations. 

• Appendix m presents the explicit expressions at leading order for the average color currents of 
partons < C >A (y 

• Appendix O com pl etes section|4]and appendix [B]by providing explicit formulae necessary to evaluate 
the MLLA corrections 5 < C >a 0 to the average color currents; 

• While the core of the paper mainly give results for LHC, Appendix 0 provides an overview at 

LEP and Tevatron energies. It is shown how, considering too large values of x (In 7 < 2) endanger 
the positivity of at low k±. Curves are also given for the range of applicability of our 

approximation is discussed in relation with the core of the paper. 

• in Appendix |E] we compare the DLA and MLLA approximations for the spectrum, the double 
differential 1 -particle inclusive distribution, and the inclusive k± distribution. 


2 THE PROCESS UNDER CONSIDERATION 


It is depicted in Fig. 1 below. In a hard collision, a parton Aq is produced, which can be a quark or a 
gluon 5 . Ho, by a succession of partonic emissions (quarks, gluons), produces a jet of opening angle ©o, 
which, in particular, contains the parton A; A splits into B and C, which hadronize respectively into the 
two hadrons h\ and /12 (and other hadrons). 0 is the angle between B and C. 

Because the virtualities of B and C are much smaller than that of A GDI. 0 can be considered to be close 
to the angle between h\ and I 12 lITOl II I II ; angular ordering is also a necessary condition for this property 
to hold. 

5 in p — p or p — p collisions, two partons collide which can create Ao either as a quark or as a gluon; in the deep inelastic 
scattering (DIS) and in e + e“ colliders, a vector boson (7 or Z) decays into a quark-antiquark pair, and Ao is a quark (or an 
antiquark); 
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Fig. 1: process under consideration: two hadrons h\ and h 2 inside one jet. 


Aq carries the energy E. With a probability Zy, it gives rise to the (virtual) parton A, which carries the 
fraction u of the energy E; (z) is the splitting function of A into B and C, carrying respectively 
the fractions uz and u(l — z) of E\ h\ carries the fraction x\ of E\ h 2 carries the fraction 1:2 of E\ 
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r . — - u( \ — z)EQ, Qo I are their respective energy distributions. 

L \w(l ~ z ) ) 

One has 0 < ©o- On the other hand, since k_ > Qo (Qo is the collinear cutoff), the emission angle 

must satisfy 0 > 0 mm = Qq/(xE), x being the fraction of the energy E earned away by this particle 

(see also subsection 12 . ll belowl. 

The following expression for the inclusive double differential 2-particle cross section has been demon¬ 
strated in m CD: 
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( 1 ) 


where 


da 


dFl 


jet / 0 


is the Born cross section for the production of Aq, Qj et is the solid angle of the jet and 


p is the azimuthal angle between B and C. 
a s (q 2 ) is the QCD running coupling constant: 

u s (q 2 ) = 


47T 


AN C pin 


q 


( 2 ) 


A 2 

1V QCD 


where A qcd ~ a few hundred MeV is the intrinsic scale of QCD and 

1 


0 = 


4 Nr 




( 3 ) 


is the first term in the perturbative expansion of the /3-function, N c is the number of colors, Tr = n//2, 
where nj is the number of light quark flavors (nj = 3); it is convenient to scale all relevant parameters 
in units of 4Aj : . 

In Q, the integrations over u and z are performed from 0 to 1; the appropriate step functions ensuring 
uz > x\, u(l — z) > X 2 (positivity of energy) are included in Dg 1 and D 
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2.1 Notations and variables 


The notations and conventions, that are used above and throughout the paper are the following. For any 
given particle with 4-momentum (ko,k), transverse momentum k± > Qq (k± is the modulus of the 
trivector k±), carrying the fraction x = k^/E of the jet energy E, one defines 


(■ = ln^ = ln(l/x), y = ln^. (4) 

ko Q o 

Qo is the infrared cutoff parameter (minimal transverse momentum). 

If the radiated parton is emitted with an angle x) with respect to the direction of the jet, one has 

k± = \k \ sin d ~ /co siri'd. (5) 


The r.h.s. of 0 uses \k\ ~ ko, resulting from the property that the virtuality k 2 of the emitted parton is 
negligible in the logarithmic approximation. For collinear emissions (t) <C 1), ~ \k\t) ~ kyd. 

One also defines the variable Y$ 


Yo = t + y = In 


E^—) 

ko QoJ 



( 6 ) 


to the opening angle @ 0 of the jet corresponds 


*0o = In 


EQ 0 

Q o ’ 


( 7 ) 


EQ o measures the “hardness” of the jet. Since ’0 < @ 0 , one has the condition, valid for any emitted soft 
parton off its “parent” 

< *0 O ■ (8) 


The partonic fragmentation function D^(xb, Q- q) represents the probability of finding the parton b hav¬ 
ing the fraction Xb of the energy of a inside the dressed parton a; the virtuality (or transverse momentum) 
k 2 a of a can go up to | Q 2 \, that of b can go down to \q 2 \. 


2.2 The jet axis 


The two quantities studied in the following paragraphs (double differential 1-particle inclusive distribu¬ 
tion and inclusive k± distribution) refer to the direction (axis) of the jet, with respect to which the angles 
are measured. We identify it with the direction of the energy flow. 

The double differential 1-particle inclusive distribution dx c ‘~ d ^ Q e is accordingly defined by summing the 
inclusive double differential 2-particle cross section over all h 2 hadrons and integrating it over their 
energy fraction x-y with a weight which is the energy (a; 2 ) itself, it measures the angular distribution of an 
outgoing hadron h\ with energy fraction x\ of the jet energy, produced at an angle 0 with respect to the 
direction of the energy flow. 

Once the axis has been fixed, a second (unweighted) integration with respect to the energy of the other 
hadron (aq) leads to the inclusive k± distribution ,[ /A , . 
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DOUBLE DIFFERENTIAL 1-PARTICLE INCLUSIVE DISTRIBU- 

d 2 N 


TION 


dx i d In © 


After integrating trivially over the azimuthal angle (at this approximation the cross-section does not 


depend on it), and going to small 0, the positive quantity dx f (Bn e 
d?N ^ f 1 , da 


reads 


dx i din 0 


V] / dX2 X2 
h 2 J ° 


d^ljet dx i dx 2 d In 0 / da 


dfl 


'jet J 0 


( 9 ) 


We use the energy conservation sum rule ITU 

V [ dxxD^(x ,...) = 1 
, Jo 


( 10 ) 


expressing that all partons /12 within a dressed parton (C) carry the total momentum of C, then make the 
change of variable v = u ^i_ z ^ where u(l — z) is the upper kinematic limit for x- 2 , to get 

V / dx 2 x 2 D } g ( X2 u(l - z)EQ,Qo] = u 2 (l - z) 2 , (11) 

^7 do \ U A~ Z ) / 

and finally obtain the desired quantity; 

d?N 


dx 1 din 0 


Y,J du J dz 1 ^ Z)D *> {U,m °' uEQ) D h B l (^,uzEe,Q 0 ) ■ 


the summation index C has been suppressed since knowing A and B fixes C. 
We can transform G3 by using the following trick: 

dz _ , f , [ dz /' ,, N r du 


du 


-(1 — z) = J du J - J d(uz ) / 
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( 12 ) 
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and o becomes 

d?N 
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£ / du Da 0 i u > EQq, uEQ) ^ f ,uzEQ,Q 0 ) 

^ d B J z TT UZ 

-Y1 [ d ( uz ) D B (?z’ uzE@ ’Qo)Y1 [ 
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u An 
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(14) 


We then make use of the two complementary DGLAP (see also the beginning of section 0 evolution 
equations fl3l which contain the Sudakov form factors dA and dp; of the partons A and B respectively: 

7 —1 /7 2 \ d I, pnn U 

d A r A ^ D A \—’ uEe ’ 


x - f dz 

1 —EJ 

B J 


$ B a (z)D h B 1 (^-,uzE@ 
\uz 


7 J 

(15) 
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the variable uz occurring in O has been introduced; in and JUT ( uEQ ) 2 refers respectively to 
the virtualities k\ and k B of A and B. Using (fl5l and (1 1 61 . (TT-il) transforms into 


= £/ duDl(u,Ee 0 ,u E e)i- A Hkl)* (^mQo) 

A ^ 

d 


(17) 


D h A depends on the virtuality of A through the variable |7] = £(k A )— £ (Qq) = 

and elementary kinematic considerations 01 lead to k A ~ (uEQ) 2 . 

By renaming B —> A and w —> u, o finally becomes 

d 2 N 


+ £/ dwD% (^,wEQ,Q 0 ) («>,£©<),«>£©)] • 

qcd)\ 


■ In 


47V c /3 l \b.(Qq/Aq CD ) J 


dx i din 0 
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F^FOq^F©^ 1 ^)-^ d A (k 2 A )D h A ' (^,uEQ,Q 0 ) 


+-°A 
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/M f£l 
u 
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din© 


din© L 

,«£©,Qo) [^(4)^1 («,F0 o ,uF0)] 

duD Ao (u, E® o, uF0) D* 1 (^, rtF0, Qq) 


( 18 ) 


and one gets 


d 2 N 


d v /ii 


( Xl ,Q,E ,® 0 ) 


with 


dxidln© din© 
i^o fall ©> E 1 ©o) = £ f duDi 0 (u,E@ 0 ,uE®)D h / (^, Uj E?0,Qo) 


(19) 

( 20 ) 


F defined in (1201 is the inclusive double differential distribution in aq and 0 with respect to the energy 
flux (the energy fraction of the hadron hi within the registered energy flux) and is represented by the 
convolution of the two functions D Aq and D A . 

The general formula (E§ is valid for all x\\ its analytical expression in the small x\ region will be written 
in the next section. 


d 2 N 

4 SOFT APPROXIMATION (SMALL-au) FOR- 

d £i dlnk± 

At t\ fixed, since y\ = ln(Z_i_/Qo) an< l Y = ln(F0/Qo) = l\ + yi, dy\ = dlnk± = din© and we 
write hereafter M f d £ k± or instead of 

Since the (/-integral {H is dominated by u = 0(1) 6 , the DGLAP jT| partonic distributions D A (u. ... ) 
are to be used and, since, on the other hand, we restrict to small x\, x\/u <C 1 and the MLLA inclusive 
D h A ((x\/u). i ...) are requested. The latter will be taken as the exact solution (see 00) of the (MLLA) 

6 D h A (^-,uEQ,Qo) « (u/xi)x (slowly varying function) - see 1221 - and the most singular possible behavior of 
Da 0 (m, EQ o, uEQ, Qo ), which could enhance the contribution of small u, is ~ 1/m; however, the integrand then behaves like 
Const, x (slowly varying function) and the contribution of small u to the integral is still negligible. 
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evolution equations that we briefly sketch out, for the sake of completeness, in appendix |A] MLLA evo¬ 
lution equations accounts for the constraints of angular ordering (like DLA but unlike DGLAP equations) 
and of energy-momentum conservation (unlike DLA). 

For soft hadrons, the behavior of the function i) h A (x\ , E(~). Qo) at x,\ <C 1 is [fl 

D h /(x 1 ,EQ,Q 0 )^-p h A 1 (ln-M^- = Y e ) , (21) 

where p A is a slowly varying function of two logarithmic variables that describes the “hump-backed” 
plateau. 

For D 1 ^ ,uEO,Q 0 ) occurring in this yields 

D h l Qo) « —p h ^ fin— ,\nu + Ye] . (22) 

V U / X~y \ X ^ / 

Because of {6}, one has 

p h A (£,Y e ) = p%(e,e + y) = !>%(£, y), (23) 

and, in what follows, we shall always consider the functions 

xD a {x, EQ,Q 0 ) = Da((-i y). (24) 


The expansion of p A ( In —,In u + Yq ) around u = 1 (In u 


u 


In 1) reads 


^-D h l ,uE@,Qo) 

such that 


Pa (4 + In u, Ye + In u ) = p \^ (4 + In u, yi + h + In u) 

Da (4 + In u, yi ) = E>a (4 ,y 1 )+\nu^-D h A 1 (h,yi) + ... ,(25) 


xiF%(x i,0,E,9 o ) ~ £ / duuDi 0 (u,EG 0 ,uE&) (b h A H,l x ,y x )+\nu dD A^ y ^ 


A L 


/ (IuuDa 0 {u, EQq,uEQ) 


Da 1 (turn) 


+ / duu\nuDA 0 (u, EQq, uEQ) 

A 


dD A (£i,yi) 

dt i 


(26) 


the second line in {26} is the 0(1 ) main contribution; the third line, which accounts for the derivatives, 
including the variation of a s , makes up corrections of relative order O(^fal) with respect to the leading 
terms (see also {37}), which have never been considered before; since, in the last line of (l26l) . u < 1 =>• 

d£) h l 

In u < 0 and d( f( is positive (see appendix I A. 41 ). the corresponding correction is negative. A detailed 
discussion of all corrections is made in subsections 14. Il and l4. 41 

It is important for further calculations that {20} has now factorized. 

While (1201) (1261) involve (inclusive) hadronic fragmentation functions D h A = I) b 1 or Dq 1 , the MLLA 
partonic functions D b A (£,y) satisfy the evolution equations (14(4 with exact solution {52}, demonstrated 
in S3 and recalled in appendix lAl The link between the latter (D q , D q , Dg, D q ) and the former goes 
as follows. At small x, since quarks are secondary products of gluons, for a given “parent”, the number 
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of emitted quarks is a universal function of the number of emitted gluons: the upper indices of emitted 
partons are thus correlated, and we can replace in (Ell the inclusive fragmentation functions by the 
partonic ones, go to the functions Da(£, y), where the upper index (which we will omit) is indifferently 
g or q, and rewrite 

e,£,e 0 ) ~ ( < u >1 +<5< u >^ Q il>AA(h,yi))DA(h,yi), (27) 

A 


with 


< U >A 0 


5<u>i 0 


f duuD^ 0 (u, EOo,uEO) « f cIuuDa 0 (u, EQq, EQ), 

J o Jo 

/ du (u \o.u)Da 0 (u, EQo, uEQ) « / du( u\o.u)Da q (u, EQq, EQ) , (28) 
Jo " Jo 


and 


tpA,h(h,yi) = 


dD A (h,yi) 


DA(h,yi) 


dE 


Thus, for a gluon jet 

xi( Xl ,Q,E,e Q ) 


~ < u > 9 g D g (£i,yi)+ <u> q g D q (h,yi) 
+ 5< u > 9 V , ff ^i(^i,yi)-Dg(^t,yi) 

+ 5<u> q g 'i/jq } t 1 (£ 1 ,yi)D q (£ 1 ,yi), 


and for a quark jet 


xiFq 1 ( Xl ,e,E,e 0 ) 


~ < u >1 Dg{h,yi)+ < U > q q D q (£i,yi) 

+ S<u>^ij}g j £ 1 (£i,yi)D g (£ 1 ,yi) 

+ 5< u > q 2/t)• 


(29) 


(30) 


(31) 


It turns out (see (TJ) that the MLLA corrections to the formulae 


Cf ~ ~ Cf ~ 

D 9 n « —D 9 D q « —D q 
q N c 9 q N c 9 ’ 


(32) 


do not modify the results and we use EJ in the following. We rewrite accordingly (Hi and (Eli 


xiF ^ (xi, 0, E, 0 O ) 
(xi,0,£;,©o) 


c* >2 tj < c ~ < c > 0 ~ 

" —^ ^ 29 g (^, yi ), 


N r 


<C>°„ +5< C >„ 


q i ” ^ ^ p, ^ _ < C > q 


N r , 


D g (h,yi) = 


N r , 


D g (h,yi), (33) 


with 

<c >° g = <U> 9 N c + <u> q C F , 

< C > q = <u> 9 N c + < u > q C F , (34) 

7 In J2SI . u is integrated form 0 to 1, while, kinematically, it cannot get lower than xi; since we are working at small xi, this 
approximation is reasonable. 








and where we have called 


5<C> g = N c S< u > 9 g ipg/^ijyi) + Cf 5< u > 9 ipq t e 1 (£i,yi), 
5<C> q = N c S< u > 9 q ^i(4,yi) + C F 5< u > q q i> q ^{h,yi)- 


< C >a 0 is the average color current of partons caught by the calorimeter. 
Plugging dTTt into dl9l ) yields the general formula 


f d 2 N \ 
\dlidln k±) q g 


d 

dyi 


<C> 


q,g 


N r 


D g (h,yi) 


(35) 


(36) 


The first line of (l30l and (ED are the leading terms, the second and third lines arc corrections. Their 
relative order is easily determined by the following relations (see d47t for the definition of 70 ) 


d?N 


1 lvr £ ^Nfuy 1 ) + j- c D a ( ei ,y,)N<c>^ 


dl 1 din A;_l 

-^-D g (£i,yi) = o(rfo) =o(y/N- s ), 

^-<c > q ,g= O( 7o 2 ) = o(a a y, 


(37) 


The different contributions are discussed in subsections 14. 1 1 and 14.41 below. 

• 7 JnTj 5 ~ = dDg J^’ y ^ (see the beginning of this section) occumng in dm is plotted in Fig. 12 and 13 of 
appendix IaI and dDa ^' y ^ occurring in (1271 ) ( 1291 ) is plotted in Figs. 14 and 15. 


• The expressions for the leading terms of x\(x\. 0, E, 0q) together with the ones of < C > { g and 
< C >° q are given in appendix iBl 


• The calculations of <5 < C > g and 5 < C > q are detailed in appendix O where the explicit analytical 
expressions for the < u >’s and d < u >’s arc also given. 


We call “naive” the approach” in which one disregards the evolution of the jet between ©o and 0; this 
amounts to taking to zero the derivative of < C > q .g in EU; EU, El, dm then yield 

< c > naive = ^ < fj > naive = Cp _ ( 38) 


4.1 The average color current < C >a 0 

On Fig. 2 below, we plot, for Y© 0 = 7.5, < C >°, < C >° q +5 < C > q , < C >°, < C >° g +5 < C > g 

as functions of y, for l = 2.5 on the left and l = 3.5 on the right. Since 0 < ©o, the curves stop at y such 
that y + ( = Y© 0 ; they reach then their respective asymptotic values N c for < C > n and Cf for < C > q , 
at which 5 < C > q and 5 < C > g also vanish (see also the naive approach El)- These corrections also 
vanish at y = 0 because they are proportional to the logarithmic derivative (1 /D(£,y))(dD(£,y)/d£) 
(see El) which both vanish, for q and 7 , at y = 0 (see appendix lAl and Figs. 16-17); there, the values 
of < C > g and < C > q can be determined from E1E1- 

The curves corresponding to LEP and Tevatron working conditions, Y© 0 = 5.2, are shown in appendix 

0 
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u o 
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Fig. 2: < C > { \ Q and < C > ( , j 1q +() < C >a 0 for quark and gluon jets, as functions of y, 
for Yq 0 = 7.5, t = 2.5 on the left and l = 3.5 on the right. 

Two types of MLLA corrections arise in our calculation, which are easily visualized on Fig. 2: 

* through the expansion (El around u = 1, the average color current < C >° 4q gets modified by 

5 < C >, 4 0 < 0 of relative order it is represented on Fig. 2 by the vertical difference between 

the straight lines (< C >^ 0 ) and the curved ones (< C >^ o +c) < C > 4 0 ); 

* the derivative of < C >° A(> with respect to y is itself of relative order Oi ^faf ) with respect to that of 
D g \ it is the slopes of the straight lines in Fig. 2. 

The y derivatives of < C >° 0 +<5 < C >a 0 differ from the ones of the leading < C > ( 4o ; this effect 
combines the two types of MLLA corrections mentioned above: the derivation of < C > with respect to 
y and the existence of 5 < C >. 

For Yq 0 = 7.5, the S < C > correction can represent 50% of < C > g at l = 2.5 and y ~ 1.5; for higher 
values of l (smaller x), as can be seen on the right figure, its importance decreases; it is remarkable that, 
when 5 < C > is large, the corrections to with respect to d< ^ ! f > become small, and vice-versa: at 

both extremities of the curves for the color current, the 5 < C > corrections vanish, but their slopes are 
very different from the ones of the straight lines corresponding to < C >°. 

So, all corrections that we have uncovered are potentially large, even dS< 2y > , which is the y derivative of 
a MLLA corrections. This raises the question of the validity of our calculations. Several conditions need 
to be fulfilled at the same time: 

* one must stay in the perturbative regime, which needs y\ > 1 (k± > 2.72A qcd ~ -7GeV; this 
condition excludes in particular the zone of very large increase of d ^^ kl _ w ^ cn 2/i —>• 0 (this property 
is linked to the divergence of the running coupling constant of QCD a s (fei) —► oo when k± —► A qcd)- 

* x must be small, that is i large enough, since this is the limit at which we have obtained analytical 
results; we see on Fig. 2 that it cannot go reasonably below l = 2.5; this lower threshold turns out to be 
of the same order magnitude as the one found in the forthcoming study of 2-particle correlations inside 
one jet in the MLLA approximation (5J; 

* (MLLA) corrections to the leading behavior must stay under control (be small “enough”); if one only 
looks at the size of the 5 < C > corrections at Y© 0 = 7.5, it would be very tempting to exclude 
y G [.5, 2.5]; however this is without taking into account the y derivatives of < C >, which also play an 
important role, as stressed above; our attitude, which will be confirmed or not by experimental results, is 
to only globally constrain the overall size of all corrections by setting x small enough. 
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Would the corrections become excessively large, the expansion (El should be pushed one step further, 
which corresponds to next-to-MLLA (NMLLA) corrections; this should then be associated with NMLLA 
evolution equations for the inclusive spectrum, which lies out of the scope of the present work. 

Though 5 < C > can be large, specially at small values of i, the positivity of < C >° +5 < C > is 
always preserved on the whole allowed range of y. 

The difference between the naive and MLLA calculations lies in neglecting or not the evolution of the 
jet between @ 0 and 0, or, in practice, in considering or not the average color current < C >a 0 as a 
constant. 

We present below our results for a gluon and for a quark jet. We choose two values Yq 0 = 7.5, which can 
be associated with the LHC environment 8 , and the unrealistic Y© 0 = 10 (see appendix [D]for Yq 0 = 5.2 
and 5.6, corresponding to the LEP and Tevatron working conditions). For each value of Yq 0 we make 
the plots for two values of ly, and compare one of them with the naive approach. 

In the rest of the paper we always consider the limiting case Q q —» A qcd 44 A ~ 0, 

A = ln-^_. (39) 

J^QCD 

The curves stop at their kinematic limit yimax such that y\ max + t\ = Yq 0 . 

d 2 N 

4.2 -at small cci: gluon jet 

diy d In k± 6 j 


On Fig. 3 below is plotted the double differential distribution of a parton inside a gluon jet as 

a function of y\ for different values of () (fixed). 


= 7.5 


; 10.0 




Fi S • 3: dS^kl foragluonjet. 


On Fig. 4 are compared, for a given value of (i, the two following cases: 

* the first corresponds to the full formulae (EHH; 

8 Sharing equally the 14 TeV of available center of mass energy between the six constituent partons of the two colliding 
nucleons yields E ~ 2.3 TeV by colliding parton, one considers a jet opening angle of 0 ss .25 and Q o « A qcd ~ 250 MeV; 
this gives Y = In « 7.7. 
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* the second corresponds to the naive approach (see the definition above El) 


d 2 N 


dl i din k _l 




( 40 ) 


dD g (h,yi) . . ■ rzzi 

—--is given in (1561) . 

dy i 




F/g. 4: 


difdMk±_ f or a S^ uon j et at fixed £\, MLLA and naive approach. 


The raise of the distribution at large k_ is due to the positive corrections already mentioned in the 
beginning of this section, which arise from the evolution of the jet between 0 and 0 q. 


d 2 N 

4.3 -at small * 1 : quark jet 

dt\ d In k± H J 


On Fig. 5 is plotted the double differential distribution of a parton inside a quark jet as a 

function of y\ for different values of t\ (fixed). 


= 7.5 


10.0 




y 


for a quark jet. 


On Fig. 6 are compared, for a given i\ fixed, the full formulae dUEl and the naive approach 


( d 2 N y aive _c F d n 

\d£idlnk±) q N c dy\ 9 1,2/1 


(41) 
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Y 0 =7.5 
u 0 


Y 0 = 10.0 
'-'o 




Fig. 6: ; n | (, r i |)' h . f or a quark jet at fixed t\, MLLA and naive approach. 


We note, like for gluon jets, at large y, a (smaller) increase of the distribution, due to taking into account 
the jet evolution between 0 and ©o- 

4.4 Comments 

The gluon distribution is always larger than the quark distribution; this can also be traced in Fig. 2 which 
measures in particular the ratio of the color currents < C > 9 / < C > q . 

The curves for d p 1 c fihik ± have been drawn for l\ = ln(l/.X] ) > 2.5; going below this threshold exposes 
to excessively large MLLA corrections. 

The signs of the two types of MLLA corrections pointed at in subsection pm vary with y: 5 < C > 
always brings a negative correction to < C >°, and to dt ' 1 ^ k ; for y > 1.5, the slope of < C > is 
always larger that the one of < C >°, while for y < 1.5 it is the opposite. It is accordingly not surprising 
that, on Figs. 4 and 6, the relative positions of the curves corresponding to the MLLA calculation and 
to a naive calculation change with the value of y. At large y, one gets a growing behavior of d ^i^ lk 
for gluon jets (Fig. 4), and a slowly decreasing one for quark jets (Fig. 6), which could not have been 
anticipated a priori. 

We study in appendix IE. 21 how MLLA results compare with DLA d D3, in which the running of a s 
has been “factored out”. 

dN 

5 INCLUSIVE k± DISTRIBUTION- 

d\nk ± 


Another quantity of interest is the inclusive k± distribution which is defined by 

( dN \ f , ( d 2 N \ _ [ Y °o-y ^ ( d 2 N \ 

\dlnkJ gorq -J Xl \dx 1 dlnkJ gorq -l min 1 {de 1 dlnk ± ) gorq ’ ^ 

it measures the transverse momentum distribution of one particle with respect to the direction of the 
energy flow (jet axis). 
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We have introduced in d42l ) a lower bound of integration i m i n because our calculations are valid for 
small x\, that is for large Aj . In a first step we take 7 rnin = 0, then vary it to study the sensitivity of the 
calculation to the region of large x\. 

We plot below the inclusive Ay distributions for gluon and quark jets, for the same two values Y@ 0 = 7.5 
and Y© 0 = 10 as above, and compare them, on the same graphs, with the “naive calculations” of the 
same quantity. 


5.1 Gluon jet; £ m in = 0 




y y 

Fig. 7: for a gluon jet, MLLA and naive approach, 

for lmm=o, Yq o = 7.5 and Y 0O = 10. 
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5.2 Quark jet; i min = 0 




Fig. 9: d ^ for a quark jet, MLLA and naive approach, 
for lmm=0, Yq o = 7.5 and Yq 0 = 10. 




Fig. 10: enlargements of Fig. 9 at large k± 


5.3 Role of the lower limit of integration £ m in 

To get an estimate of the sensitivity of the calculation of d ^ k to the lower bound of integration in l42l) . 
we plot in Fig. 11 below the two results obtained at Yq 0 = 7.5 for t min = 2 and i mtn = 0, for a gluon 
jet (left) and a quark jet (right). 
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Y 0 =7.5 
u 0 


Y e =7.5 
'-'0 




Fig- lit t j ] T j j. with (’min — 2 and (min — 0 
for gluon (left) and quark (right) jet. 


The shapes of the corresponding distributions are identical; they only differ by a vertical shift which is 
small in the perturbative region y > 1 (restricting the domain of integration - increasing ( rnm - results 
as expected in a decrease of ). This shows that, though our calculation is only valid at small x \, 
the sensitivity of the final result to this parameter is small. 

5.4 Discussion 

MLLA corrections are seen on Fig. 8 and Fig. 10 to cure the problems of positivity which occur in the 
naive approach. 

The range of (\ integration in the definition (El) of d should be such that, at least, its upper bound 
corresponds to x \ small enough; we have seen in the discussion of MLLA corrections to the color current 
in subsection Olthat one should reasonably consider t\ > 2.5; at fixed Yq 0 this yields the upper bound 
2/t < >©o — 2.5, that is, at LHC y\ < 5. 

On the other side, the perturbative regime we suppose to start at y\ > 1. These mark the limits of 
the interval where our calculation can be trusted 1 < y\ < 5 at LHC. For y i < 1 non-perturbative 
corrections will dominate, and for y\ > Y© 0 — (""" ~ Yq 0 — 2.5, the integration defining rl ranges 
over values of x\ which lie outside our small x approximation and for which the MLLA corrections 
become accordingly out of control. 

On the curves of Figs. 7 and 9 at Y(- )(} = 10, the small y region exhibits a bump which comes from the 
competition between two phenomena: the divergence of a s (k^_) when k± —> Q o and coherence effects 
which deplete multiple production at very small momentum. The separation of these two effects is still 
more visible at Yq q = 15, which is studied in appendix IE.3I where a comparison with DLA calculations 
is performed. At smaller Y© 0 , the divergence of a s wins over coherence effects and the bump disappears. 
The curves corresponding to the LEP and Tevatron working conditions are given in appendix lol 

5.4.1 Mixed quark and gluon jets 


In many experiments, the nature of the jet (quark or gluon) is not determined, and one simply detects 
outgoing hadrons, which can originate from either type; one then introduces a “mixing” parameter u. 
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which is to be determined experimentally, such that, for example if one deals with the inclusive k± 
distribution 

/ dN \ ( dN \ . . ( dN \ 

V d In k± ) mixed UJ \dlnkx) g + U \dlnk ± ) q ' 

It is in this framework that forthcoming data from the LHC will be compared with our theoretical pre¬ 
dictions; since outgoing charged hadrons are detected, one introduces the phenomenological parameter 
JC ch (2 (7| normalizing partonic distributions to the ones of charged hadrons 


dN r rA dN ) 

d In k± J \ d In k± ) 

/ \ / mixed 


(44) 


6 CONCLUSION 


After deducing a general formula, valid for all x, for the double differential 2-particle inclusive cross 
section for jet production in a hard collision process, the exact solutions of the MLLA evolution equa¬ 
tions 0 have been used to perform a small x calculation of the double differential 1-particle inclusive 
distributions and of the inclusive k± distributions for quark and gluon jets. 

Sizable differences with the naive approach in which one forgets the jet evolution between its opening 
angle ©o and the emission angle 0 have been found; their role is emphasized to recover, in particular - , 
the positivity of the distributions. 

MLLA corrections increase with x and decrease when the transverse momentum k± of the outgoing 
hadrons gets larger; that they stay “within control” requires in practice that the small x region should 
not be extended beyond £ < 2.5; it is remarkable that similar bounds arise in the study of 2-particle 
correlations 0. At fixed Y© 0 , the lower bound for i translates into an upper bound for y: this fixes in 
particular the upper limit of confidence for our calculation of q j^' ; above this threshold, though k± is 
larger (more “perturbative”), the small x approximation is no longer valid. 

The “divergent” behavior of the MLLA distributions for y —► 0 forbids extending the confidence domain 
of MLLA lower that y > 1, keeping away from the singularity of a s (k[ j_) when k± — > Aqcd- 

The two (competing) effects of coherence (damping of multiple production at small momentum) and 
divergence of a s (k^_) at small k± for the inclusive k± distribution have been exhibited. 

MLLA and DLA calculations have been compared; in “modified” MLLA calculations, we have further¬ 
more factored out the a s dependence to ease the comparison with DLA. 

While the goal of this work is a comparison of our theoretical predictions with forthcoming data from 
LHC and Tevatron, we have also given results for LEP. LHC energies will provide a larger trustable 
domain of comparison with theoretical predictions at small x. 

Further developments of this work aim at getting rid of the limit Qo ~ A qcd and extending the calcu¬ 
lations to a larger range of values of x; then, because of the lack of analytical expressions, the general 
formulae <m and <EB should be numerically investigated, which will also provide a deeper insight into 
the connection between DGLAP and MLLA evolution equations m. 

Acknowledgments: It is a pleasure to thank M. Cacciari, Yu.L. Dokshitzer and G.P. Salam for many 
stimulating discussions, and for expert help in numerical calculations. R. P-R. wants to specially thank 
Y.L. Dokshitzer for his guidance and encouragements. 
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APPENDIX 


A EXACT SOLUTION OF THE MLLA EVOLUTION EQUATION FOR 
THE FRAGMENTATION FUNCTIONS; THE SPECTRUM AND ITS 
DERIVATIVES 

A.l MLLA evolution equation for a gluon jet 


Because of IS, we will only write the evolution equations for gluonic fragmentation functions D b . 

The partonic structure functions D b satisfy an evolution equation which is best written when expressed 
in terms of the variables l and y and the functions D b defined by (TJ (see also (ED<E1): 

XbD b a (x b , k a , q) = D b a (£ b , y b ). (45) 


The parton content D g of a gluon is shown in 00 to satisfy the evolution equation (Y and y are linked by 

©) 


D, 


,(£, y) = 5(£) + [ V dy' [ + y') [l - a5(t’-£)\ D g {£', y'), (46) 

Jo Jo 


'0 j 0 

where the anomalous dimension 70 (y) is given by (A is defined in d39l ) 


7o (V) = 4A C 


a 


:( k l) 


1 


2 vr f3(y + X) 

(see the beginning of section[ 2 ]for 13, Tr, Cf, a s , N c ) and 


a = 


AN r 


11 


-N c + -T R 1 - 


2 C F 

IvT 


; C F = A/3 for SU(3) C . 


(47) 


(48) 


The (single logarithmic) subtraction term proportional to a in © accounts for gluon —> quark transi¬ 
tions in parton cascades as well as for energy conservation - the so-called “hard corrections” to parton 
cascading 

No superscript has been written in the structure functions D g because the same equation is valid indif¬ 
ferently for Dg and D g (see section 0). One considers that the same evolution equations govern the 
(inclusive) hadronic distributions D 1 / (Local Hadron Parton Duality). 


A.2 Exact solution of the MLLA evolution equation for particle spectra 


The exact solution of the evolution equation ©, which includes constraints of energy conservation and 
the running of a s , is demonstrated in |Sj to be given by the following Mellin’s representation 


D g (£,y, A) 





ds 

v + s 


\(u + s)vj 


d ,l/ a )t+vy 

2n i 


V + S J 


a/0 


e 


—As 
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(49) 











From d49l) and taking the high energy limit £ + y = Y A 9 one gets (I] |(7j the explicit formula 


Dg(t,v) = 


+ V f %e~^^{A + l,B + 2M^ + y)), 


( 50 ) 


/ 3B(B + 1) J 27rz 

where $ is the confluent hypergeometric function the integral representation of which reads G3 d 


$(A + 1,5 + 2, w7) = T(5 + 2) (wF) 


-s-t 


dt t 


-B 


with -4 = + B = f 


J (27ri) — 1) — 1 

/*oo 

r (n) = / dxx n_ 1 e _x . 

■Jo 


uYt. 


(51) 


Exchanging the t and cu integrations of <HB (ED and going from t to the new variable a = ^ In 
El becomes 


Dg(t, y ) = ^ 2 ^ + B (r, y, l) j , 

1 y 

where the integration is performed with respect to r defined by a = - In — + it, 

2 £ 


^B{r,y,£) = 


cosh a — - -- sinh a 

_ y + t 

£ + y a 


1 B/2 


f3 sinh a 


I B (2y/Z(T,y,t)), 


2 t - 1 
(52) 


Z(r,y,£) = ~~~j~ ■ t fcosh a — - sinh a ) ; 

p smh a \ y + £ 


Ib is the modified Bessel function of the first kind. 


(53) 


A.3 The spectrum 


On Fig. 12 below, we represent, on the left, the spectrum as a function of the transverse momentum (via 
y) for fixed i and, on the right, as a function of the energy (via £) for fixed transverse momentum. 




Fig. 12: spectrum D(£,y) of emitted partons 
as functions of transverse momentum (left) and energy (right) 
9 Y EQ Qo/Fqcd is not strictly equivalent to Q o —> A qcd (limiting spectrum). 
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Fig. 13 shows enlargements of Fig. 12 for small values of y and £ respectively; they ease the understand¬ 
ing of the curves for the derivatives of the spectrum presented in subsection IA.4I 




Fig. 13: spectrum D(t. y) of emitted partons 
as functions of transverse momentum (left) and energy (right): enlargement of Fig. 12 


A comparison between MLLA and DLA calculations of the spectrum is done in appendix IE. II 


A.4 Derivatives of the spectrum 


We evaluate below the derivatives of the spectrum w.r.t. In k± and ln( 1 /x). 

We make use of the following property for the confluent hypergeometric functions <f> fiTH : 

— (A + 1, B + 2, u (£ + y)) = — (A + 1, B + 2, u> (£ + y)) = to 4> (A + 2, B + 3, u (l + y )) . 
di dy B + 2 

(54) 

• We first determine the derivative w.r.t. £ = ln(l jx). Differentiating d50l w.r.t. £, and expanding (El. 
one gets 10 |H 


d_ 

dl 


D g (£, y) = 2 


T(B) /•! dr 

P Jo ^ 


—Ba 


—J— (1 + 2e“ sinh a) T B + \ e °^B+\ 
£ ~r y fj 


(55) 


k± 

Differentiating w.r.t. y = In — yields 

Q o 


d 


dy 9 ~ (3 


D g {l,y) = 2 


T(B) ["2 dr —Ba. 


7T 


jJ— (1 + 2e“ sinha) T B + \e a T B +i -^^F B -i 
£+ y P l+y 


( 56 ) 


In Fig. 14, Fig. 15, Fig. 16 and Fig. 17 below, we draw the curves for: 

dD g (£, y) 


dy 

dD g (£,y ) 

dy 


as a function of y, for different values of l fixed; 
as a function of l, for different values of y fixed; 


" <551 and <56> have also been checked by numerically differentiating <52t . 
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* — g f„ ’ ^ as a function of £ for different values of y fixed; 


di 

dDg{£, y ) 
d£ 


as a function of y for different values of l fixed. 


In each case the right figure is an enlargement, close to the origin of axes, of the left figure. 




Fig. 14: 


dDg(C,y) 

dy 


as a function of y for different values ofi 




Fig. 15: 


dD g (t,y) 

dy 


as a function of £ for different values ofy 




Fig. 16: 


dP B (t,V) 

dt 


as a function of £ for different values ofy 
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Fig. 17: 


dDg(£,y) 

clt 


as a function of y for different values of£ 


That 


dD g (l,y) 

dy 


derivative. 


goes to infinity when y —> 0 is in agreement with the analytic behavior in In (£/y) of this 


B LEADING CONTRIBUTIONS TO x x ( x u &, E , 0 O ) AT SMALL 

Xi 


Using (f32t . the leading terms of x i (x\ , 0, E, 0q) (l26l calculated at small x± read 

xiFg 1 (xi,0,U,0 o )° ~ Dg(£i,yi) (< u > 9 g < u = < N >9 -P g (li,yi), 

XlFq 1 (xi,0,U,0 o )° « Dg(£i,yi) (<u>°+^-<u>{\ 

(57) 

The leading < C >° and < C >° q in (134b for a quark and a gluon jet are given respectively by (see fll . 
chapt. 9 11 ) 


< C> 


o 

9 


< C >oo —ci (N c — Cf) 

< C >oo —ci ( N c — Cf) 


( In (EQ/ Aqcd) \ (c3W) 
V ln (EQ 0 / Aqcd) ) 

/ y© + a \ (cs/4JVc/3) 

U>„ + V 


(58) 


<C> 


o 

a 


< C >oo +C2 (Nc — Cf) 


< C >oo +C2 ( N c — Cp) 


/ In (EQ/ Aqcd) \ (c3/4jVe/?) 
\ln (EQ 0 / Aqcd) ) 

/ y© + a \ (cs/4iVc/3) 

U>„ + V 


(59) 


u The coefficient /3, omitted in the exponents of eqs. (9.12a), (9.12b), (9.12c) of [TJ has been restored here. The factor 4iV c 
is due to our normalization (see the beginning of section^. 
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with 


<C> oo = CiN c + C2Cf, 

Cl 

in the r.h.s of m 
oo, 0 —> 0. 

In practice, we take in this work 

Qo ~ A qcd o A ~ 0, (61) 

which ensures in particular the consistency with the analytical calculation of the MLLA spectrum (ap¬ 
pendix [At. which can only be explicitly achieved in this limit. 


8 Cf 
3^"’ 


1 2 nr % 2 

c 2 = 1 c\ — — c 3 = -C F + -nf, 
3 c 3 3 3 


(60) 


we have used the definitions @ 0. < C > oo corresponds to the limit E 


C CALCULATION OF 8 < C > g and S < C > q OF SECTION! 

C.l Explicit expressions for < u >^ o and S < u >^ o defined in (l28l) 


The expressions d28l) for < u and 5<u >(J (j are conveniently obtained from the Mellin-transformed 
DGLAP fragmentation functions 0 


V(j,0 = [ duu J 1 D(u,£), 

Jo 

which, if one deals with D^(u, r 2 , s 2 ), depends in reality on the difference £(r 2 ) — £(s 2 ): 


m 2 ) = f 

J Li 


Q 2 dk 2 a s (k 2 ) 
k 2 47 t 


e(r 2 )-^ 2 ) 


■ In 


' \n(r 2 / JJqc D ) 


4 N c (3 \ ln(s 2 /Al CD ) ) ' 


One has accordingly 


(62) 


(63) 


<u>i = vi 0 (2,a(E& 0 ) - am), s< u a = -vi 0 (j,aE @ 0 ) - am) 


3 =2 


(64) 


The DGLAP functions V(j, £) are expressed 0 in terms of the anomalous dimensions vrij), i' gU) an d 
u±(j), the j dependence of which is in particular known. 

For the sake of completeness, we give below the expressions for the < u >’s and d < u >’s. 


< u > 


g 

a 


JL /7 ke 0 + A \i \ ( Y eo + A \~i 
25 y\Y e + \J J \y© + a y 5 


< u > 


9 

9 


1/25 



( Ye o + A A 
\Ye + \) 



( Y&o + A 
V A e + A 


50 

81 


< u > 9 q 


16 // y eo + A \i \ ( Ye 0 + A \-t 
25 [\Ye + \J J \ 1© + A / ’ 
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< u 


< u 


< u > s q ea + < u 

5 < 


5 < 


6 < 


5 < u 


> 

n 


= -1/25 -9 


Yeo± AV 1 
+© + A 


-16 + 25 


^0o+A \ 2/9N 

Yq + XJ 


Yq q + A 
le + A 


50 

81 


> 


val 


+0 O + A 
+© + A 


> 


val 


= 1/25 9 


+0 O + A 
le + A 


+ 16 


+00 + A 
le + A 


u>l 


1 


337500 


-43011 


+6804 7T 2 - 48600 In 
+0o + A 


+21600 In 


-17400 In 


+© + A 

+©o + A 

+0 + A 


) 50 
81 

+0o + A 
+0 + A 
+00 + A 


+ 43011 - 6804 vr 2 

V 50 

+0o+ A ^ 

+0 + A 


50 

81 


7 n 


+e + A 

+0q + A 
+0 + A 


vr 2 + 109525 In 


_ 50 
81 


+0q+A ^BI 

+0 + A 


+©o + A 

+© + A 


u>9 9 


1 


337500 


-86400 In 


-11664 

+0q + A 
,+0 + A 

+0o + A 


+0q + A 
+© + A 


+ 31104 vr 2 


V 50 

+0O+A V 1 
le + A 
+©o + A 


+38400 b ^ . 2 

V+0 + A/ \ +© + A / ^ 

+11664 - 31104 vr 2 - 109525 In | 100 + A 

V +© + A 

'+0O + A\ 2 \ { +©o + A 


+17400 In 


+© + A 


7r 


+© + A 


+e 0 + A A 81 
le + A 


u > 9 q 


759375 


48114 


+0o + A yi 
+© + A 


-48114- 6804 7T 2 


x 50 

+00 + AA 81 


+© + A 


+6804 vr 2 - 48600 In 


+21600 In 


-17400 In 


+0q + A 
+0 + A 
+0q + A 
+0 + A 


+©o Y X\ ( +0 O + A ^ 81 
+© + A 


+© + A 
+00 + A 


50 

81 


7 r 


+© + A 
+0q + A 
+© + A 


vr 2 + 109525 In 


_50 

81 


+ 0 Q + A 
+0 + A 


> 
^ n 


759375 


-13122 


100 + 81 + 34992 vr 2 

+© + A 


+0q + A 
+© + A 
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+54675 In 


-24300 In 


*e„ + A 
F© + A 
Yq 0 + A 


F© + A / \ Y e + A 


50 

^eo + A yi 
F© + A 
>00 + A 


7r 


+13122 - 34992 vr 2 + 219050 In 


-265625 In 


+37500 In 


F+ + A 
. F© + A 
>e 0 + A 


>b 0 + A 
F© + A 
>0o + A A 2/9 
Y e + Xj 
10 O + A\ 2/9 


Yq + X ) \ 7e + A / 


- 34800 In 



^0o + A 

F© + A 


Yq 0 + A 
F© + A 


7r 


50 

81 


5 < u > 


= -4j(- 85 + 12 ’' 2 ) ln (tTX 


> 0p + A 

F© + A 


jr ^ ^ val i c- ^ ^ sea 

0 < U > +0 < U > q 


759375 


-54675 In 


+24300 In 


-219050 In 


13122 


V 50 

*e 0 + A yi 
F© + A 


- 34992 


7T 


2 ( F© 0 + A 
F© + A 


>0o + A A / F© 0 + A 


Je + A 
^0o + A 
, F© + A 
{ Yq q + A 
Ue + A 


F© + A 
Yq 0 + A 
. Y e + A 

+ 34800 In 


50 
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When © —> ©o, all <5 < it >’s vanish, ensuring that the limits £(.E©o) — £{EQ) —> 0 of the « c >° Ao 
+<5 < C >a 0 )’ s are the same as the ones of the < C >^ Q ’s. 


C.2 8 < C > q and 8< C > g 


They are given in and one uses d32l such that only ip g .e } (see d29t ) appears. Their full analytical 
expressions for the 8 < C >’s are too complicated to be easily written and manipulated. 

Using the formulae of IC.ll one gets the approximate results 


6<C> c 


1.4676 - 1.4676 


+0.5461 



(4,l/i), (66) 


and 
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<5< C > c 


-2.1898 + 2.1898 [ Y ®° ^ - 3.2510 In 


-0.3072 


Yq + X 

T© 0 + X 
Yq + A 


.50 s 

In ( Y&0 + A ) | 'ipg^ii^yi). 


Y e + A 


( 67 ) 


The logarithmic derivative '(/++ ( + 1 , J/ i ) d29l) of the MLLA spectrum D g [l\ , r/i) is obtained from ll52l) of 
appendix lAl 


D AT LEP AND TEYATRON 


At LEP energy, the working conditions correspond to Yq 0 ~ 5.2; they are not very different at the 
Tevatron where T© 0 ~ 5.6. We first present the curves for LEP, then, after the discussion concerning 
the size of the corrections and the domain of validity of our calculations, we give our predictions for the 
Tevatron. 


D.l The average color current 




Fig. 18 < C >^ o and < C >^ Q +<5 < C >A 0 for quark and gluon jets , as functions of y, 
for Tq 0 = 5.2, £ = 1.5 on the left and l = 2.5 on the right. 


Owing to the size of the (MLLA) corrections to the < C >’s and their y derivatives, we will keep to the 
lower bound l\ > 2.5. 


D.2 


d 2 N 

dlx d In k± 


for a gluon jet 


We plot below for the two values 


1.5 and t = 2.5. 
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Y a . = 5.2 


Y„. = 5.2 




Fig. 19: for a gluon jet affixed t \, MLLA and naive approach. 

The excessive size of the 5< C > corrections emphasized in subsection [Dj] translates here into the loss 
of the positivity for d ^ d ^ lk± at £ = 1.5 for y < 1: our approximation is clearly not trustable there. 

d 2 N 


D.3 


for a quark jet 


dt i d In k 

We consider the same two values of t as above. 


Y ra . = 5.2 


Y„. = 5.2 




Fig. 20: , 1 (^, 1 ^ f.. fi >r a quark jet affixed l\, MLLA and naive approach. 
Like for the gluon jet, we encounter positivity problems at l = 1.5 for y < 1.25. 


dN 

D.4 -for a gluon jet 

d In k i 


We plot below for a gluon jet obtained by the “naive” approach and including the jet evolution 

from 0o to 0; on the right is an enlargement which shows how positivity is recovered when MLLA 
corrections are included. 
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Fig. 21: for a gluon jet, MLLA and naive approach. 


dN 

D.5 -for a quark jet 

dlnfcj^ 


We proceed like for a gluon jet. The curves below show the restoration of positivity by MLLA correc¬ 
tions. 


Y„ = 5.2 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

0 1 2 3 4 5 

y 

Fig. 22: for a quark jet, MLLA and naive approach. 

That the upper bound of the £\ domain of integration defining d ^ corresponds to a large enough 
t\ > 2.5 requires that, for LEP, y\ should be smaller that 5.2 — 2.5 = 2.7; combined with the necessity 
to stay in the perturbative regime, it yields 1 < Vi < 2.7. 

D.6 Discussion and predictions for the Tevatron 

The similar condition at Tevatron is 1 < y\ < 5.6 — 2.5 = 3.1; like for LEP, it does not extend to large 
values of k± because, there, the small x approximation is no longer valid. We give below the curves that 
we predict in this confidence interval. 



2.5 3 3.5 4 4.5 5 

y 
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Y 0 =5.6 
u 0 


Y 0 =5.6 




Fig. 23: d i^l ± for a gluon (left) and a quark (right) jets, MLLA predictions for the Tevatron. 

Since experimental results involve a mixture of gluon and quark jets, the mixing parameter lo (subsection 
EZD has to be introduced in the comparison with theoretical curves, together with the phenomenological 
factor K ch normalizing partonic to charge hadrons distributions. 


E COMPARING DLA AND MLLA APPROXIMATIONS 


DLA El El and MLLA approximations are very different Q); in particular, the exact balance of 
energy (recoil effects of partons) is not accounted for in DLA. 

We compare DLA and MLLA results for the two distributions of concern in this work. Studying first 
their difference for the spectrum itself eases the rest of the comparison. 

We choose the two values Yq 0 = 7.5 and Yq 0 = 15. While the first corresponds to the LHC working 
conditions (see footnote |HJ», the second is purely academic since, taking for example ©o ~ .5 and Qq ~ 
250 MeV, it corresponds to an energy of 1635 TeV ; it is however suitable, as we shall see in subsection 
IE.3l to disentangle the effects of coherence and the ones of the divergence of a s at low energy in the 
calculation of the inclusive k± distribution. 

E.l The spectrum 


Fixing a s in DLA at the largest scale of the process, the collision energy, enormously damps the cor¬ 
responding spectrum (it does not take into account the growing of a s accompanying parton cascading), 
which gives an unrealistic aspect to the comparison. 

This is why, as far as the spectra are concerned, we shall compare their MLLA evaluation with that 
obtained from the latter by taking to zero the coefficient a given in <HT which also entails B = 0; 
Fo(T,yJ) in <|53} becomes /o(2\/ Z(r, y.F). The infinite normalization that occurs in d52l because of 
F(B = 0) we replace by a constant such that the two calculations can be easily compared. This realizes 
a DLA approximation (no accounting for recoil effects) “with running a s ”. 

On Fig. 24 below are plotted the spectrum D g (£, y = Y© 0 — £) for gluon jets in the MLLA and DLA 
“with running a s ” approximations. 
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Fig. 24: the spectrum D 9 (£,Yq 0 — t) for gluon jets; 
comparison between MLLA and DLA (“with running a s ”) calculations. 


The peak of the MLLA spectrum is seen, as expected, to occur at smaller values of the energy than that 
of DLA. 


E.2 Double differential 1-particle inclusive distribution 


The genuine MLLA calculations being already shown on Figs. 3 and 5, Fig. 25 displays, on the left, 
a “modified” MLLA calculation obtained by dividing by a s (k\ ) ~ 2 n py ( see with A —> 0); 
subtracting in the MLLA calculations the dependence on k± due to the running of a s (k‘j_) allows a 
better comparison with DLA (with fixed a s ) by getting rid of the divergence when k± —> Qq. 

On the right are plotted the DLA results for gluon jets, in which a s has been fixed at the collision energy 
(it is thus very small). Since their normalizations are now different, only the shapes of the two types of 
curves must be compared; we indeed observe that the DLA growing of with k± (or y \) also 

occurs in the “modified” MLLA curves. 




Fig. 25: comparison between MLLA (after dividing by a:. s (k j ), on the left) 
and DLA calculation with a s fixed (on the right) of d ^ f or gluon jets. 


The DLA distribution for quark jets is obtained from that of gluon jets by multiplication by the factor 
Cf/N c ; it it thus also a growing function of y\. 
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The MLLA distribution for quark jets, which is, unlike that for gluon jets, a decreasing function of y\ 
(see Fig. 6), becomes, like the latter, growing, after the dependence on a s (k‘j_) has been factored out: 
one finds the same behavior as in DLA. 


E.3 Inclusive k± distribution 


On Fig. 26 we have plotted, at Yq 0 = 7.5: 

- the MLLA calculation of d ^ divided by a s (k‘j_), such that the divergence due to the running of a s 
has been factored out, leaving unperturbed the damping due to coherence effects; 

- the DLA calculation of d with a s fixed at the collision energy. 

Like in lE.2l because of the division by a s , the two curves are not normalized alike, such that only their 
shapes should be compared. 

The comparison of the DLA curve (at fixed a s ) with the genuine MLLA calculation displayed in Fig. 7 
(left) shows how different are the outputs of the two approximations; while at large k± they are both 
decreasing, at small k± the running of a s makes the sole MLLA distribution diverge when k± —► Q o 
(non-perturbative domain). 




y y 

Fig. 26: Y eo = 7.5: comparing MLLA and DLA calculations °f -jTnFf ( see a ^ so ^8- 7); 
from left to right: — f^-MLLA and DLA (a s fixed). 

C*S\rCj') 


In the extremely high domain of energy Yq 0 = 15 used for Fig. 27, the two competing phenomena 
occurring at small y\ can then be neatly distinguished. 

The first plot, showing MLLA results, cleanly separates coherence effects from the running of a s \ in the 
second figure we have plotted the MLLA calculation divided by a s (k“j_): damping at small y\ due to 
coherence effects appears now unspoiled; finally, DLA calculations clearly exhibit, too, the damping due 
to coherence 12 . 

The large difference of magnitude observed between the first (genuine MLLA) and the last (DLA) plots 
occurs because DLA calculations have been performed with a s fixed at the very high collision energy. 

Like in lE~2l because of the division by a s , the second curve is not normalized like the two others, such 
that only its shape should be compared with theirs. 

12 The DLA points corresponding to y\ = 0 can be analytically determined to be 4 N c /nf (gluon jet) and 4 Cf/ rif (quark 
jet); they are independent of the energy Te 0 . 
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Fig. 27: Y &0 = 15: comparing MLLA and DLA calculations of d ^ ; 

from left to right: MLLA , — jr^-MLLA and DLA (a s fixed). 

ot s yfcj,) 
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Figure captions 


Fig. 1: the process under consideration: two hadrons h\ and h -2 inside one jet; 

Fig. 2: < C >^ Q and < C >^ Q +5 < C >a 0 for quark and gluon jets, as functions of y, for Y© 0 = 7.5, 
t = 2.5 and £ = 3.5; 

Fig. 3: dtfd Efcl for a g luon jet, Yq 0 = 7.5 and Y © 0 = 10; 

Fig. 4: d(:i d j^ lk± at fixed i\ for a gluon jet, comparison between MLLA and the naive approach; 

Fig. 5: dh d dhi k ± for a quark jet, Y© 0 = 7.5 and Y 0O = 10; 

Fig. 6: d £ d d ^ lk± at fixed £\ for a quark jet, comparison between MLLA and the naive approach; 

Fig. 7: inclusive k± distribution d ^ for a gluon jet, Y© 0 = 7.5 and Y© 0 = 10, MLLA and naive 
approach, both for £ nan = 0; 

Fig. 8: enlargements of Fig. 6 at large k±; 

Fig. 9: inclusive k± distribution d ^ for a quark jet, Y© 0 = 7.5 and Y© 0 = 10, MLLA and naive 
approach, both for £ nan = 0; 

Fig. 10: enlargements of Fig. 8 at large k±\ 

Fig. 11: role of the upper limit of integration over x\ in the inclusive k± distribution d ^ for gluon 
(left) and quark (right) jet; 

Fig. 12: spectrum D g (£,y) of emitted pardons as functions of transverse momentum (left) and energy 
(right); 

Fig. 13: enlargements of Fig. 11 close to the origin; 

Fig. 14: d ^ 9 d y' y ' 1 as a function of y for different values of t, 

Fig. 15: d ^ a }y' v ^ as a function of 7 for different values of y\ 

Fig. 16: as a f unc tion of i for different values of y\ 

Fig. 17: d ^aYy) as a f unc tion of y for different values of 

Fig. 18: < C >° 4() and < C >° t| . +5 < C > ,\ (j for quark and gluon jets, as functions of y, for Yq 0 = 5.2, 
t = 1.5 and £ = 2.5; 

Fig. 19: for a gluon jet for Y© 0 = 5.2 at fixed i\, MLLA and naive approach; 

Fig. 20: for a quar'k jet for Yq 0 = 5.2 at fixed i\, MLLA and naive approach; 

Fig. 21: d d ^, for a gluon jet for Y© 0 = 5.2, MLLA and naive approach; 

Fig. 22: d for a quar'k jet for 1© 0 = 5.2, MLLA and naive approach; 

Fig. 23: d d ^ k± for a gluon and a quar'k jets, MLLA predictions for the Tevatron. 

Fig. 24: the spectrum D g (£, Y 0O —l) for gluon jets; comparison between MLLA and DLA (“with running 
«/’) calculations; 

Fig. 25: comparison between MLLA (after dividing by a s {k! 2 1 ), on the left) and DLA calculation with 
a s fixed (on the right) of dy d ^ lk ^ for gluon jets; 
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Fig. 26: y© 0 = 7.5: comparing MLLA and DLA calculations of d ^ (see also Fig. 6); from left to 
right: —-W-MLLA and DLA (a s fixed); 

Fig. 27: Yq 0 = 15: comparing MLLA and DLA calculations of from left to right: MLLA, 

—tWtMLLA and DLA (a s fixed). 

&s\kj') 
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